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Abstract 

The luminosity of fading type la supernovae is governed by radioactive decays of ^®Ni and ^®Co. 
The decay rates are proportional to the Fermi coupling constant Gp and, therefore, are determined 
by the vacuum expectation value v of the Brout-Englert-Higgs field. We use publicly available sets 
of lightcurves of type la supernova at various redshifts to constrain possible spacetime variations of 
the ®®Ni decay rate. The resulting constraint is not very tight; however, it is the only direct bound 
on the variation of the decay rate for redshifts up to z ~ 1. We discuss potential applications of 
the result to searches for non-constancy of Gp and v. 
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I. INTRODUCTION. 


Probing the limits of applicability of fundamental laws of Nature is one of noble tasks for a 
present-day physicist or astrophysicist. In particular, tests of constancy of parameters of the 
Standard Model of particle physics and of General Relativity, dubbed fundamental physical 
constants, attract interest for decades (see e.g. Refs. BE] for recent reviews and extensive 
lists of references). Astrophysics and cosmology, working with largest distances and longest 
time periods available for research, provide for excellent opportunities to constrain spacetime 
variations of the fundamental constants. 

One of the tools to perform these studies is to compare particular astrophysical processes 
taking place at various distances and in various directions. Here, we concentrate on the 
latest stages of explosions of type la supernovae. Being very bright, these exploding stars 
have been observed at cosmologically large distances. At the same time, studies of relatively 
nearby explosions revealed clear universality in features of these supernovae and related 
their observable parameters to underlying physical processes. While precise details of the 
explosion itself remain to be understood, the fading luminosity of a supernova days to 
months after the explosion is governed by beta decays of ®®Ni and ^®Co, as it has been 
conhrmed in 2014 by direct observations 011]. This fact has been used 0 0 to constrain 
spacetime variations of the decay rates, which are determined by the Fermi constant Gp. In 
the Standard Model, this “constant” is directly related to the vacuum expectation value v 
of the Brout-Englert-Higgs scalar held. This value is determined by the effective potential 
for the scalar and may, in principle, vary with the evolution of the Universe. 

In this work, we take advantage of the fact that type la supernovae have been extensively 
studied in recent years because of their importance as distance indicators in cosmology. 
We use the most recent JLA publicly available set 0, see http://supernovae.in2p3.fr/ 
sdss_snls_jla/ReadMe .html, The light curves are available for 740 objects from various 
observational campaigns, spanning redshifts between 0 and ~ 1.3. This enlarges signihcantly 
both the number and the redshift range of supernovae, as compared to previous studies. This 
allows us to obtain the hrst constraints on possible spacetime variations of the ^®Ni decay 
rate for the redshift range up to 2 ; ~ 1. They may be related to variations in Gp and v under 
certain assumptions. 

The relation between the amount of decaying isotopes and the supernova luminosity is 
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not direct and we refer the reader to previous studies mu for discussions of this point. 
The bottom line is that after the maximum, the luminosity fading is well described by two 
exponential functions, related to the ^®Ni and ^®Co decays, with the slopes A proportional to 
the decay rates T with a coefficient of order one. Given considerable universality in all type 
la supernova explosions, the coefficient of proportionality does not vary signihcantly from 
one object to another and constraints on the variations in the brightness slope are directly 
translated into constraints on the variations of T. 


II. ANALYSIS AND RESULTS. 

We are dealing here with high-redshift objects, so the observations, performed with a 
certain hlter, correspond to various emitted wavelengths. The relation between the luminos¬ 
ity and the amount of the isotope may be wavelength-dependent, thus causing a systematic 
redshift-dependent bias in our study. For each supernova, we considered the light curve 
measured in the R hlter and calculated the corresponding emitted wavelength, given known 
z. The light curve was then htted by a standard wavelength-dependent template obtained 
by interpolation of templates of Ref. |H] available at https: //c3. Ibl. gov/nugent/nugent_ 
templates .html, Example supernova light curves are presented in Fig. together with our 
£t. 

It is well known that type la supernova light curves are not precisely universal; there 
exists a scatter in the time dependence of the hux parametrized by the so-called stretch 
factor s, one of the parameters of the £t. It is probably determined by variations in the 
initial chemical composition of the exploding star. This stretch factor relates directly to the 
quantity we are interested here, the decay rate, since it rescales the time axis. Therefore, the 
slope of the decay exponent may be determined by approximating a segment of the standard 
light curve by an exponential function, Rst('?‘) cx: exp(—AstT), where F is the hux and r is 
the time after the maximum. We choose the B-hlter template and select the hrst 30 days 
after maximum to obtain Ast — 0.0964 day“^. Then, for a supernova with a stretch factor s, 
we have A = sAst (note that the trivial redshift multiplicative correction, (1 -|- z), is already 
included at the stage of htting). 

For our analysis, we selected supernovae with at least hve data points in the light curve. 
We removed objects for which the ht was poor. There remained 687 objects in the sample 
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FIG. 1. Example light curves in the R band: SNLS-04Dlhy {z = 0.85; top panel) and SN 1995ac 
{z = 0.049; bottom panel). The lines represent best hts obtained using redshift-corrected templates. 

after applying these cuts. Figure [^presents the redshift-corrected measured lightcurve slopes 
A versus the redshift. 

In this hgure, one can see the intrinsic scatter of the decay slopes, reflecting variations in 
chemical composition, absorption etc. It is manifested by the fact that the data points at a 
given redshift are scattered far beyond their error bars. We account for this intrinsic scatter 
by introducing a universal systematic error, added in quadrature to the statistical error of 
the slope measurement. The value of the systematic error is adjusted in such a way that the 
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FIG. 2. The lightcurve slope A versus the redshift z for the sample of supernovae considered in 
this work. 

best-fit of the A( 2 ;) dependence takes the median value expected for the corresponding 
number of degrees of freedoirQ Assuming linear variation, r( 2 ;) = r(0) -|- const x x, we 
obtain the best-fit value of the slope in the redshift range 0 < z < 1.0, 

r(2; ~ 1.0) — r(x = 0) 


r(z = 0) 


= -0.182 ±0.021. 


We will return to the interpretation of this would-be signihcant effect in Sec. |III| here we 
derive a conservative constraint from it, 

r( 2 ; ~ 1.0) — T{z = 0) 


< 0.223 (95% CL). 


(1) 


r(j = 0) 

Next, we make use of the fact that the observed objects are located in different parts of 
the sky and search for global patterns of directional variations of L. To this end, we study 
the distribution of the lightcurve slopes A on the sky and fit it with a dipole harmonics in 
the equatorial coordinates (a, 5), 

X = Xq + ux n(a, 5) ■ n(ao, %), 


^ Note that for the relevant number of degrees of freedom, there is almost no difference between median 
and average. 


5 





















































Isotope 

Method 

Result 

Redshift Reference 

56Ni 

supernovae 

> -0.223 

1.0 

this work 

5®Co 

supernovae 

-0.272 ±0.156 

0.024 

El 


meteorites 

-0.016 ±0.016 

0.45 

m 


TABLE I. Comparison of astrophysical constraints on variations of nuclear beta decay rates. 


where n(a, is a unit vector pointing to the direction with coordinates {a, h); the parameters 
of the £t are the monopole harmonic Aq, the direction n(ao, 5o) and the magnitude a\ of the 
dipole. The dipole anisotropy is constrained as 

^ < 0.148 (95% CL). (2) 

Ao 

Alternatively, the first harmonic modulation in the right ascension a, 

A = Ao + cos (a — ao), 

is constrained as 

^ < 0.051 (95% CL). (3) 

Ao 

The constraints ([^, (|^, (|^ are the main direct results of our study. 


III. DISCUSSION AND CONCLUSIONS. 

The constraints on the variation of the ®®Ni decay rate obtained here may be compared to 
other (scarce) limits on the variation of beta decay rates, see Tablej^ The next step, that is to 
relate the observed result to the constancy of fundamental constants, is not straightforward 
and requires additional assumptions (see e.g. an excellent description of how contrived this 
step is in Ref. [10], Sec. 2.2, not to mention numerous more general discussions in the 
literature). 

First, the variations may be constrained in a self-consistent way for dimensionless con¬ 
stants only, while the decay rates and related, more fundamental, parameters Cp and v are 
all dimensionful. A plausible dimensionless combination for our analysis may be formed 
with the help of the gravitational constant C, cf. the “/xh” constant of Ref. m- This raises 
the second subtle point, which constants are allowed to vary. For instance, allowing for 
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variations of the gravitational constant may change the conclusions signihcantly. Therefore, 
constraining just the variation of a constant is not that useful; a particular physical model 
predicting the variations of the entire set of fundamental constants should be constrained in¬ 
stead. The results presented here, that is constraints on the variations of a physical quantity, 
the decay rate, may be used in any work of that kind. 

However, just to put our results in context and to compare them with other limits, we 
interpret our bounds on the variation of T in terms of v within a simplihed assumption that 
no other constants are allowed to vary. This is presented for demonstration only, because 
we do not mean any particular underlying theory behind the assumption. 

Apart from the approach we use here, there exist other ways to constrain variations of v. 
The change of another radioactive decay rate, also proportional to Gp, was constrained 
with meteorite data. People analyzed the effect of v on the Big Bang nucleosynthesis 
(BBN, through several quantities) and on the Cosmic Microwave Background (CMB) and 
the Large-Scale Structure (LSS) of the Universe (through the change in the electron mass 
which is proportional to n). All indirect constraints require model-dependent interpretation 
which is especially hard when effects of the strong interaction, which are not calculable in 
the fundamental theory, are taken into account. For completeness, we mention here also 
the possibility to constrain variations of v through the change in the electron mass rUe by 
constraining variations of the electron-proton mass ratio /i, the most restrictive constraints 
coming from measurements with atomic clocks (AC) and with quasar (QSO) spectra. These 
result may constrain variations of v if it is assumed that the proton mass does not change 
with the change of v while the electorn mass does, because the former is dominated by 
strong-coupling effects and is not given by the sum of quark masses, which are in turn 
proportional to v in the high-energy limit. This point requires, however, a quantitative 
derivation which is presently missing. A far from being exhaustive list of recent constraints 
and references is presented in Table [IT] 

Studies of the light-curve shapes of type la supernovae, in particular of their redshift de¬ 
pendence, have a long story. Probably the hrst work on this subject was one by Pskovskii [22] . 
Studies of the rise time of type la supernovae, whose redshift dependence was discussed, in 
particular, in Ref. [23], can hardly be directly used for the purposes of the present work. 
Modern studies of the fading part of the light curve, e.g. 1231 . suggest a trend to decreasing 
of the decay slope of the curve as the redshift increases, similarly to what we observe in the 
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Quantity 

Method 

Result 

Redshift Reference 

Av/v 

supernovae 

< 0.056 

1.0 

this work 


supernovae 

0.068 ± 0.039 

0.024 

m 


meteorites 

0.004 ± 0.004 

0.45 

m 


BBN 

0.0134 ±0.0014 

~ 10® 

m 


BBN 

0.004 ± 0.002 

~ 10® 

m 

Arue/me 

CMB 

0.004 ±0.011 

~ 10® 

m 


LSS 

U.UUD_g Q^g 

~ 10® 

mi 

A/i//i 

QSO 

(0.0 ±1.0) X 10-'^ 

0.89 

mi 


QSO 

(-2.4 ± 2.3) X 10"^ 

0.89 

m 


QSO 

(0.0 ±1.5) X 10-® 

1.3 

m 


QSO 

(-7.6 ±8.1) X 10-6 

2.4 

m 


QSO 

(-1.7 ± 1.7) X 10-6 

3.17 

m 

A/h 

AC 

(-0.5 ± 1.6) X 10-16 yj.-i 

0 

m 


AC 

(0.2 ±1.1) X 10-16 yr-i 

0 

m 


TABLE II. Comparison of recent constraints on variations of v, nie and obtained by various 
methods; see text for details. 

present work. A possible explanation is that the slope depends on the environment in the 
supernova host galaxy: the decay of the light curve is slower in young stellar environments, 
while the age of stellar populations decreases with increasing redshift. Therefore, selecting 
only supernovae residing in galaxies with active star formation may result in a more homo¬ 
geneous sample [25] and thus reduce the corresponding bias in the study of the decay rates. 
This is a prospective way to improve the precision of constraints presented here. 

To conclude, we have obtained constraints on the spacetime variation of the ®®Ni decay 
rate; our results are given in Eqs. ([^, (|^, (|^. These are the first bounds on the variation 
of the nuclear decay rates obtained for the redshift range as large as 0 < 2 : < 1, see Table 
for comparison with other results. The interpretation of the results in terms of varying 
fundamental constants like the Higgs scalar vacuum expectation value v may require a 
formulation of a particular theory describing how these, and other, constants vary. Within 










simplified assumptions, we compared our result in terms of the variation of v to others, 
including more model-dependent ones, in Table [ITI 
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